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ABSTRACT 


A crucial feature not widely accounted for in local helioseismology is that surface 
magnetic regions actually open a window from the interior into the solar atmosphere, 
and that the seismic waves leak through this window, reflect high in the atmosphere, 
and then re-enter the interior to rejoin the seismic wave field normally confined there. 
In a series of recent numerical studies using translation invariant atmospheres, we 
utilised a “directional time-distance helioseismology” measurement scheme to study 
the implications of the returning fast and Alf ven waves higher up in the solar atmo¬ 
spher e on the seismology at the photosphere ( [Cally fc Moradi||2013 Moradi & Cally 
20141. In this study, we extend our directional time-distance analysis to more realis- 


tic sunspot-like atmospheres to better understand the direct effects of the magnetic 
field on helioseismic travel-time measurements in sunspots. In line with our previ¬ 
ous findings, we uncover a distinct frequency-dependant directional behaviour in the 
travel-time measurements, consistent with the signatures of MHD mode conversion. 
We found this to be the case regardless of the sunspot field strength or depth of its 
Wilson depression. We also isolated and analysed the direct contribution from purely 
thermal perturbations to the measured travel times, finding that waves propagating in 
the umbra are much more sensitive to the underlying thermal effects of the sunspot. 


Key words: Sun: helioseismology - Sun: oscillations - Sun: surface magnetism 


1 INTRODUCTION 

Sunspots and active regions (magnetic flux concentrations 
tens of thousands of kilometres across containing sunspots) 
are the most visible manifestation of solar magnetic activity 
on the solar surface. A detailed understanding of sunspots 
and magnetically active regions is therefore essential in order 
to establish accurate physical relationships between internal 
solar properties and magnetic activity in the photosphere. 

Using observations of surface oscillations, helioseismol¬ 
ogy provides the most effective way to observationally probe 
structure inside the Sun. The combination of high spatial 
resolution, continuous observing, and simultaneous vector 
magnetograms provided by the Helioseismic and Magnetic 
Imager (HMI) instrument on board the Solar Dynamics Ob¬ 
servatory (SDO) delivers unprecedented probing of magnetic 
active regions and sunspots. 

However, important developments in the techniques of 
local helioseismology (i.e., both theoretical and observa¬ 
tional) are required to realize the full potential that these 
observations offer. As highlighted by a number of detailed 
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comparative studies and reviews on the matter (e.g., Gi- 
zon et al. 2009 Moradi et al. 2010 Moradi 20121, major 


challenges exist in the development of new helioseismic pro¬ 
cedures that are robust in the presence of magnetism and 
capable of probing both subsurface magnetic structures and 
associated flows. 

Over the years various local helioseismic techniques 
have substantially contributed to our understanding of the 


solar interior (see Gizon, Birch & Spruit 2010 for a com¬ 
prehensive review). The most widely used measurement 
method in local helioseismology is time-distance helioseis- 
mology ( Duvall et al.|[l9^ h By correlating observations of 
Doppler velocity at different times and positions on the so¬ 
lar surface, a causative link is inferred and a travel time 
between pairs of points putatively determined. Comparing 
these travel times with those calculated for the quiet-Sun, 
one infers the presence of wave-speed anomalies beneath the 
surface that may be due to such features as magnetic helds, 
temperature variations or plasma flows. 

With the noise level being substantially larger for group¬ 


time measurements (Kosovichev, Duvall & Scherrer 20001, 


time-distance travel times are typically derived from phase 
travel-times. While this is adequate for the quiet Sun, in- 
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terpretation becomes more complicated when considering 
active regions, where phase shifts can naturally arise from 
changes in wave propagation speed (e.g., due to subsurface 
flows and sound speed perturbations induced by the pres¬ 
ence of the magnetic field), but they can also result from 


other sources as well. For example, mode damping (Woodard 


1997), acoustic source suppression (Gizon & Birch 2002), 
and the Wilson depression (Briiggen fc Spruit|2000 Lindsey 


fc Braun||2000 1 have all been identified as possible sources 


of phase shifts in active regions. 

Another important source of phase shifts is via MHD 
mode conversion in the atmosphere. Mode conversion takes 
place in regions where the sound (cs) and Alfven speed (ca) 
are comparable. It is expected to be significant for sunspot 
seismology because in the umbrae of sunspots, the layer 
where Ca = Ca lies is just a few hundred kilometres below the 
formation height of the Fe I spectral line of the Helioseismic 
and Magnetic Imager (HMI) instrument on board the Solar 
Dynamics Observatory (SDO) . 

“Fast- to-slow” mode conversion has been explored as 
a primary cause of acoustic wave (p-mode) absorption in 
sunspots for decades (jSpruit & Bogdan 1992 Cally, Bog-| 
dan fc Zweibel||1994[]CaIly||1995| [Crouch fc Cally||2003[ 


Shelyag et al.|2009 1. In this scenario, the p—modes emerging 


in sunspots below the Ca = Cs level are effectively (acoustic) 
fast waves. On passing through the Ca = Cs layer the p-modes 
are partially transmitted into the solar atmosphere as (pri¬ 
marily acoustic) slow waves, most efficiently at small “attack 
angle” (the angle between the wavevector and the magnetic 
field). The transmitted sound waves propagate longitudi¬ 
nally along field lines at frequencies above the field-adjusted 


acoustic cutoff frequency and are reflected otherwise | Bel & 
Leroy|1977 1. If the attack angle is not small however, signif¬ 


icant amounts of energy will be converted to (magnetic) fast 
waves ( Schunker fc Cail^|2006 ). These fast waves are then 
reflected off the Alfven wave speed gradient, at the height 
where their horizontal phase speed (v^h = ^Ikh', where cu is 
the angular frequency and kh the horizontal wavenumber) is 
approximately equal to Ca, back down to the surface (having 
assumed Ca S> Ca at this level). 

More recently it has been realised that fast waves cre¬ 
ated in this way are further subject to partial conversion to 
Alfven waves higher in the atmosphere. Depending on the 
local relative inclinations and orientations of the background 
magnetic field and the wavevector, the fast wave may un¬ 
dergo partial mode conversion to either an upward or down¬ 
ward propagating Alfven wave around the reflection height, 
where they are near-resonant (Cally fc Goossens|2008 Cally 


& Hansen 2011 Khomenko & Cally 2011 2012 Felipe 2012). 


After they reflect off the Alfven wave speed gradient, the 
fast waves may re-enter the solar interior wave field. This 
could be problematic for helioseismology, since any phase 
changes produced by the “fast-to-Alfven” mode conversion 
process would seriously compromise any inferences derived 


from helioseismic 

inversions of phase travel times (e.g., Du- 

vail et al.|[1996[ [Kosovichev, Duvall & Scherrer[[2000||Cou- 

vidat et al.[[2005 

1 , which would normally, but inaccurately, 


interpret such phase changes as “travel-time shifts” due to 
subsurface inhomogeneities alone. 


In a series of recent numerical studies (Cally & Moradi 
[Moradi fc CaII^[M^ [2^, we quantified the im¬ 
plications of the returning fast and Alfven waves for the 


seismology of the photosphere by comparing Alfvenic losses 
higher up in the solar atmosphere with helioseismic travel¬ 
time shifts at the surface. Using 3-D numerical simulations 
of helioseismic wave propagation in simple translationally in¬ 
variant atmospheres, we applied a “directional time-distance 
helioseismology” approach sensitive to magnetic field orien¬ 
tation, finding substantial wave “travel time” discrepancies 
of several tens of seconds (depending on field strength, fre¬ 
quency, and wavenumber) related to phase changes resulting 
from mode conversion, and not “actual” travel time changes. 
These results, which were also verified using the Boundary 


Value Problem (BVP) method of Cally & Goossens (2008) 
and Cally (2009), indicated that processes occurring higher 


up in the atmosphere are strongly influencing the core data 
products of helioseismology. 

In these studies only translation invariant setups were 
used, which are most useful to study the effect in funda¬ 
mental terms, but are not typically found on the Sun. Our 
best chance at constraining the interior structure of sunspots 
comes with constructing accurate forward models. Hence, in 
order to be able to make meaningful estimates of the direct 
role played by MHD mode conversion and wave reflection 
in helioseismic measurements, we extend our study to more 
realistic sunspot model atmospheres spanning the subphoto¬ 
sphere {z = —10 Mm, with z = 0 being the photosphere) to 
the chromosphere {z = 1.9 Mm) and study the sensitivity 
of directional helioseismology measurements to changes in 
the photospheric and subsurface structure of sunspot mod¬ 
els. For practical computational reasons, we were unable to 
model the seismic effects of the Transition Region, though 
[Hansen fc Callyj (|2014|) find it too has significant signatures. 


2 THE BACKGROUND MODEL 

The background models we employ consist of a number of 
azimuthally symmetric magneto-hydrostatic (MHS) sunspot 


atmospheres adopted from Khomenko & Collados (2008). 


To summarise, these models consist of a concatenation of a 
self-similar model in the deep photospheric layers, calculated 


following the method of Low (1980), with a potential solu¬ 


tion above some arbitrary height using the method of|P izzo[ 
(1986). The thermodynamic variables for the sunspot axis 


are taken from the semi-empirical Avrett (19811 umbral 


core model, while the “quiet-Sun” atmosphere variables are 
taken from Model S ( Christensen-Dalsgaard et al.[[19^ ) in 
the deep sub-photospheric layers, smoothly joined to the 
VAL-C ( Vernazza, Avrett fc Loeser|1981 1 model in the pho¬ 
tospheric and chromospheric layers, and stabilised using the 


method outlined by Parchevsky & Kosovichev (2007). 


The sunspot models possess a high degree of flexibility 
for conducting a detailed directional helioseismology study, 
with a number of variable parameters such as field strength 
on the axis, field inclination at the photosphere, and spot 
radius. Another important variable parameter in the mod¬ 
els, is the Wilson depression - the height difference between 
the umbra and photosphere. This can be easily changed in 
the model by choosing the desired location of the constant 
optical depth log rsooo = 0 (the formation height of the 5000 
A continuum radiation) of the semi-empirical umbral model. 
Studies have shown that the Wilson depression may be a sig- 


nifiant source of travel time reductions in sunspots (Braun 
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Figure 1. An example of the MHS sunspot atmosphere used in 
the forward modelling. A cut through y = 0 for a 1.5 kG sunspot 
with a Wilson depression of 400 km is shown here. The back¬ 
ground colour-scale corresponds to logCa (km s~^). The contour 
lines represents the field strength (in G). The dotted line indi¬ 
cates the location of the plasma ^ Pi 1 layer, while the dashed 
line represents the reference observation height corresponding to 
logT-sooo = -1-6. 


fc Lindsey|20Q0 Lindsey, Gaily fc Rempel|2010 1 , so it is im¬ 
portant to observe what effects it may have on directional 
travel-time measurements. 

We conduct a number of experiments with the sunspot 
models. In the first set of experiments, we study the sensitiv¬ 
ity of directional travel times to sunspot field strength and 
inclination using two sunspots models with differing peak 
photospheric field strengths (1.5 and 2.5 kG) but with all 
other parameters fixed (i.e., spot size, field inclination at 
the photosphere and Wilson depression, which is fixed at 
400 km). In the second set of experiments, we investigate 
the sensitivity of the directional travel times to the depth 
of the Wilson depression. Estimates from observations put 
the depth of Wilson depression in the range of 300 — 1500 


km ( 

Bray & Loughhead 

1964 Martinez Billet & Vazquez 

[1993 

[Mathew et al.|2004 

Watson et al.|2009[). For our study 

we employ three identica 

. sunspot models with a relatively 


moderate surface field strength (1.5 kG) with varying Wil¬ 
son depressions depths of 300, 400 and 500 km respectively. 
In our final set of experiments, we ascertain the contribu¬ 
tion of the underlying thermal perturbations to the travel¬ 
time shifts in contrast to the direct magnetic effects in a 
sunspot with surface field strength of 1.5 kG and a Wilson 


depression of 400 km. As shown in Moradi & Gaily (20081 
and [Moradi, Hanasoge fc Cally| ( |2009| this can easily be 
achieved in linear numerical simulations by suppressing the 
direct magnetic effect on the waves. The combined outcomes 
from these experiments provide us with valuable diagnostics 
of both the thermal and magnetic structures of sunspots. 


As all sunspot models encompass a Wilson depression, 
the simulated data are analysed at a constant optical depth 
logrsooo = — 1 . 6 , which roughly represents the layer where 
the contribution function for Fel 6173 A photospheric spec¬ 
tral line has its maximum (Khomenko et al.|2009 1. The opti¬ 
cal depths are calculated using the routine described injjessj 


et al. (20121 by integrating the continuum and line opaci¬ 


ties along the lines-of-sight for each column in the sunspot 
models. The ATLAS9 package | Kurucz||1993 1 opacities are 
used in the computation. We also have the ability to choose 
the line-of-sight viewing angle (from the vertical), but for 
simplicity we calculate and compare directional travel times 
using the photospheric velocity at disk centre. Some proper¬ 
ties of these sunspot model atmospheres are shown in Figure 

m 


3 THE FORWARD MODEL 


As in our previous studies, we numerically solve the lin¬ 
earised equations of ideal MHD using the Seismic Propaga¬ 
tion through Active Regions and Convection (SPARC) code 
( Hanasoge|2007 1 which has been successfully utilised in the 
past to study wave propagation through model sunspots 
( [Moradi, Hanasoge fc Cally||200^ . The dimensions of the 
3-D computational box employed for the numerical simu¬ 
lations are 140 Mm in the horizontal (i, y) directions and 
11.9 Mm in the vertical (z) direction. The bottom bound¬ 
ary of the domain is located at 10 Mm below the photo¬ 
spheric level 2 = 0. The horizontal grid spacing consists of 
256 equidistant points in x and y, with a resulting resolu¬ 
tion of Ax = Ay « 0.55 km/pixel, while the vertical grid 
spacing A 2 is nonuniform, ranging from tens of kilometres 
near and above the surface to just over one hundred kilome¬ 
tres near the bottom of the computational domain. The top 
~ 500 and bottom ~ 800 km of the box are occupied by the 
vertical absorbing (PML) boundary layers, while absorbing 
sponges line the sides of the box. 

The axis of the sunspot is placed at the centre of the 
computational domain. For each sunspot case studied, we 
conduct ten unique simulations using a Gaussian perturba¬ 
tion source positioned along the left hand side of the sunspot 
along y = 0, starting from the axis (x = 0 , 3 / = 0 ) and then 
at nine other locations along the negative x axis, as depicted 
in Figure As the sunspot model is axisymmetric, this al¬ 
lows us to study each corresponding 9 associated with the 
source location separately. 

The acoustic source employed for our calculations is 


similar to that employed by Shelyag et al. (20091 and Moradi 


& Gaily (20141, where a source term of the form: 


. 27rt / (r - ro) 

Vz = sm exp--— 

ti V '^r 


exp|--A_|A)^ (1) 


is added to the right hand side of the vertical momentum 
equation. In the equation above Vz is the perturbation to 
the vertical component of the velocity, to = 300s, G = 300s, 
at = 100s, Gr = 4Ax, and ro(x, y) is the source position. The 
source, which is always initiated below the surface at 2 = 
—0.65 Mm, generates a broad spectrum of acoustic waves in 
the 3.33 mHz range, mimicking wave excitation in the Sun. 
For each magnetic case we conduct a separate quiet-Sun run 
to act as the reference (unperturbed) model. 
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The relatively large field strengths associated with the 
sunspot models being considered, coupled with the exponen¬ 
tial drop in density with height in the atmosphere, naturally 
results in a substantial Ca above the surface. For explicit nu¬ 
merical solvers, such as SPARC, this results in severe CFL 
(At ~ Az/ca) constraints, significantly compounding the 
computational expense of conducting a detailed parametric 
study. To alleviate the problem, we employ a Lorentz Force 
“limiter” to limit/cap the Alfven wave speed at a particular 
value above the surface. This approach is commonly adopted 
by explicit numerical solvers in computational MHD stud¬ 
ies of sunspot structure I Rempel, Schiissler fc Knolker|20d9 


Cameron et al.|2011 Braun et al.|2012 1 , allowing one to in¬ 

crease the simulation At to any desired or practical value. 


However in |Moradi fc Cally] ( |2014[ ), we studied the phys¬ 
ical implications of imposing an artificial limit on Ca and 
found that it can severely impact the fast-wave reflection 
height (ca ~ ijJikh) in the sunspot atmosphere, which can 
be problematic for fast-to-Alfven mode conversion and any 
subsequent helioseismic analyses. In fact, we found that un¬ 
less the Ca cap is placed well above the horizontal phase 
speed associated with the wave travel distance being stud¬ 
ied (thus ensuring minimal damage to the fast-wave reflec¬ 
tion height), helioseismic travel time measurements could 
be severely affected. On the back of these findings, we de¬ 
cided to employ a limiter with a Ca cap at 80 km s“^, but 
restrict our helioseismic analyses to waves with horizontal 
phase speeds well below this (see Table [^, so as to ensure 
our travel time measurements would not be compromised. 


4 DIRECTIONAL TIME-DISTANCE 

HELIOSEISMOLOGY 

With single source wave excitation, time-distance diagrams 
can easily be constructed by plotting the resulting velocity 
signal as functions of time for all horizontal locations. More¬ 
over, each source location along the negative x-axis corre¬ 
sponds to a specific field inclination 9 (from the vertical). As 
seen in Figure with the source locations we have chosen 
we can sample 9 in the range 0° — 70°. By selecting a receiver 
location at a horizontal distance (A) away from the source 
location around the xj/-plane, we isolate the magnetic field 
orientation with respect to the vertical plane of wave prop¬ 
agation, which we refer to as the “azimuthal” field angle {(j>, 
where 0 ^ (j) ^ 180°, from the right- to the left-hand side of 
the sunspot) which we sample in 10° bins. 

Prior to calculating the travel times, we first filter the 
data cubes in two frequency ranges: 3 and 5 mHz by em¬ 
ploying a Gaussian frequency filter with a dispersion of 0.5 
mHz. We also apply an /-mode filter to remove the con¬ 
tribution from surface gravity waves. We then measure the 
phase travel time perturbations St (i.e., the differences in the 
phase travel times between the magnetic and nonmagnetic 


simulations) using Gabor wavelet fits | Kosovichev & Duvall 


19971 to the time-distance diagram at various A away from 


the source, for each source (9) receiver {(p) pair of points. 
A rectangular window of width 14 minutes ceutred on the 
first-bounce ridge selects the fitting interval in time lag. The 
fits are done by minimising the misfit between the Gabor 
wavelet and the wave form. An initial guess of the Gabor 
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Figure 2. The one-way travel time measurement geometry. In 
both panels, the background represents the vertical component 
of magnetic field strength at the observation height logrsooo = 
— 1.6. In panel a) the contours are indicative of the magnetic field 
inclination from the vertical {0 in degrees) at the same height 
in the atmosphere, while the crosses represent the locations of 
the individual acoustic sources utilised in the forward modelling 
calculations. Panel b) shows an example of the receiver locations 
(which span from 0^0^ 180°, from the right- to the left-hand 
side of the sunspot, spaced 0 = 10° apart) for a source initiated 
at 0 = 0°, denoted by the cross on the axis (x = 0, y = 0). The 
dots indicate the receiver locations at A = 6.2 Mm, the circles 
A = 8.7 Mm, and the diamonds represent A = 11.6 Mm. 


wavelet parameter values is obtained by fitting the refer¬ 
ence (quiet-Sun) wave form first. We use MATLAB’s mul¬ 
tidimensional unconstrained nonlinear minimisation routine 
fminsearch for the fitting, which employs the Nelder-Mead 
simplex algorithm ( [Lagarias, et al.||1998[ ). This is a direct 
search method that does not use numerical or analytic gra¬ 
dients. We measured 5t for three typical skip distances A 
( [Couvidat et al.||2005[ ). The horizontal phase speeds associ¬ 
ated with these distances are presented in Table 
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Table 1. Wave travel distances analysed and their associated 
horizontal phase speeds in the quiet Sun model. 


A (Mm) 

Vph (km s 1) 

6.2 

12.6 

8.7 

14.1 

11.6 

16.4 


5 RESULTS & ANALYSIS 


5.1 Sensitivity of Directional Travel Times to 
Frequency, Field Strength and Inclination 


The contour plots in Figures[^and|^depict the time-distance 
phase travel-time perturbations (with respect to the quiet 
solar model) as functions of wave source position/held in¬ 
clination from vertical 9 and receiver location/azimuthal di¬ 
rection (j), derived from sunspot models with surface held 
strengths of 1.5 and 2.5 kG respectively. The results shown 
are for the two frequency bands analysed, 3 (left column) 
and 5 (right column) mHz, and for waves which travel a 
horizontal distance of A = 6.2 (panels a-b), 8.7 (panels c-d) 
and 11.6 Mm (panels e-f) from the source. 

Pleasingly, the features we observe in Figures - f^ 
are very much in accord with the directional travel times 
derived from previous studies using simple translationally 


invariant background atmospheres (Gaily & Moradi 2013 
Moradi fc Cally||2014 1. Specihcally, they show a clear man¬ 


ifestation of the acoustic cutoff aX 9 = 30° — 40° for 5 mHz 
and 9 = 50° — 60° at 3 mHz. For 9 below the acoustic cut¬ 
off, small positive 5t values of a few seconds are apparent. 
For larger 9 (i.e, sufficient for the ramp effect to take hold 
u) > ujcCos9), the atmosphere is open to wave penetration 
and mode conversion. This results in signihcant negative 5t 
for these 9, particularly at small sin 0 (around 0° and 180°), 
which is due to the fast magnetically-dominated waves un¬ 
dergoing significant phase enhancement on returning to the 
surface after passing upward through the Ca = Cs layer, re¬ 
flecting near to/kh = Ca, and finally re-entering the interior 
via Ca = Cs again. However, away from 0 = 0 (and 180°), 
the fast waves lose energy as they are partially converted to 
the Alfven wave, which results in a phase retardation that 
partially cancels the underlying negative travel time per¬ 
turbation at small sin0. In line with previous studies, the 
energy loss is at its maximum around 0 = 80° — 100°, orien¬ 
tations typically associated with peak fast-to-Alfven conver¬ 


sion ( Khomenko fc Gally|20lf 2012 Gaily fc Moradi|2013 1. 


However unlike our previous studies where the atmo¬ 
sphere and magnetic field were horizontally invariant, the 
presence of the sunspot, coupled with the distribution of 
the individual wave sources (on the left-hand side of the 
sunspot), results in a distinct asymmetry in 5t about 0. 
This is essentially due to one end of the wave path being in 
a stronger region of perturbation (i.e., inside the sunspot), 
and the other end being near or inside the “quiet Sun” re¬ 
gion, which will naturally result in a directional bias in 5t, 
with larger (negative) St expected for 0 < 50° (i.e., waves 
travelling primarily to the right/inside the “umbra” of our 
sunspot model). This effect is exacerbated as the wave travel 
distance is increased, as is evident in Figures e) and 

f), for A = 11.6 Mm. 






0 80 160 
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Figure 3. One-way phase travel-time perturbations (<5 t) derived 
from the 1.5 kG sunspot model calculations as a function of field 
inclination (0) from the vertical, and azimuthal angle (0) for wave 
travel distances of A = 6.2 (a-b), 8.7 (c-d) and 11.6 (e-f) Mm. 
Left column represents 3 mHz and right column 5 mHz. 


As expected, the magnitude of the Sr perturbations is 
also strongly dependant on frequency and magnetic field 
strength. Larger negative and smaller positive St are ob¬ 
served for all sunspot models and A as the frequency is 
increased from 3 to 5 mHz. This frequency dependance of 
helioseismic travel times has been well documented in the 


past 

(Braun & Lindsey 

2000 

Ghou 

20001 

Braun & Birch 

2006 

ICouvidat & Rajaguru|20071|Moradi, Hanasoge & Gaily 

2009 

. Increasing the held strength of the sunspot naturally 


shifts the location of the Ca = Cs layer deeper below the 
surface, but the only direct effect on the directional St we 
observe at the photosphere is an increase in their magnitude 
at both 3 and 5 mHz for all A. 


5.2 The Effect of the Wilson Depression 

Figurej^shows the directional St derived for A = 6.2 Mm for 
three 1.5 kG sunspot models with varying Wilson depression 
depths (300, 400 and 500 km), calculated at 3 and 5 mHz. 
While the general behaviour of St across 9 and 0 for all 
three models is consistent with those derived in section [STT] 
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Figure 4. Same as Figure]^ but for the 2.5 kG sunspot model. 


it is also apparent that modifying the depth of the Wilson 
depression can have a direct and measurable impact on the 
directional travel times. 

This is not entirely unexpected of course, as the Wil¬ 
son depression is a physical displacement in the photosphere 
which will naturally give rise to a change in the path length 
of the waves. Modifying the depth of the Wilson depression 
also implies a change in the near surface density and tem¬ 
perature stratihcation of the sunspot, which in turn will also 
modify the actual wave speed (both Cs and Ca), the result of 
which should manifest itself in the travel time calculations. 

These effects are evident in Figure which shows a cut 
at 0 « 32° through Figure Here we can clearly observe 
faster travel times associated with waves travelling towards 
the sunspot axis as the Wilson depression is shifted deeper 
below the surface. While the St differences between 300 and 
400 km are very subtle (under ~ 1 second), more significant 
differences in St are observed with the Wilson depression 
at 500 km. At <)) = 0 we see a ~ 3 second difference to the 
300 — 400 km cases, and at 5 mHz it’s ~ 7 seconds. Waves 
travelling away from spot centre (large (j>) do not appear to 
be affected by the change in Wilson depression depth. We 
observed a similar behaviour for A = 8.7 and 11.7 Mm (but 
have not shown them here for the sake of brevity). 


Figure 5. One-way phase travel time perturbations St as a func¬ 
tion of field inclination from the vertical 8 azimuthal angle <j> for 
A = 6.2 Mm, derived from three 1.5 kG sunspot models with 
Wilson depression of 300 (a-b), 400 (c-d) and 500 (e-f) km. Left 
column represents 3 mHz and right column 5 mHz. 


These results are generally consistent with the recent 
findings of Schunker et al. (20131, who studied the sen¬ 
sitivity of helioseismic travel times to the depth of the 
Wilson depression using numerical forward modelling of 
plane wave packets through non-MHS sunspot model atmo¬ 
spheres. They found that a ~ 50 km change in the Wilson 
depression can be detected above the observational noise 
level. 


5.3 The Effect of Thermal Perturbations 


One of the advantages of forward modelling of waves in a 
model sunspot atmosphere is that it provides us with the 
opportunity to isolate the individual effects of the mag¬ 
netic field and thermal perturbations on travel-time mea¬ 


surements (Moradi & Cally 2008 Moradi, Hanasoge & Cally 
|2009p . In order to isolate the thermal contributions to the 
measured directional St, we repeat our single source calcu¬ 
lations, this time using a “thermal” sunspot model, where 
only the thermal perturbations corresponding to the 1.5 
kG sunspot model with a Wilson depression of 400 km are 
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Figure 6. One-way phase travel time perturbations St for waves 
initiated at 6 ^ 32° as a function of azimuthal angle for A = 
6.2 Mm, derived from three 1.5 kG sunspot models with varying 
Wilson depression depths. The solid lines are the results from 300 
km model, the dotted lines represents the 400 km model, and the 
dashed lines represent the 500 km model. Panel a) represents 3 
mHz travel times and panel b) 5 mHz travel times. 


present, but with the direct magnetic effects on the waves 
suppressed. The directional travel times are then measured 
in an identical manner as before, with the results shown in 
Figure These travel times can be directly compared with 
those derived from the 1.5 kG magnetic sunspot in Figure]^ 
which we do so in Figure where we show some line plots 
with the thermal and magnetic travel times plotted on the 
same scale for a selection of 9. 

It is important to note that the resulting “thermal 
travel-time perturbations” produced from these calculations 
result from a combination of thermal perturbations and ge¬ 
ometrical effects due to the presence of a Wilson depression. 
When considering purely thermal effects on their own, i.e, 
a cooler plasma with a reduced sound speed and no Wil¬ 
son depression, one would expect to see positive travel-time 
shifts with respect to the quiet Sun, as waves travel slower 
in the cooler medium. On the other hand, the presence of 
a Wilson depression can change the wave-path length, de¬ 
pending on the wave propagation direction (towards or way 
from the umbra for example) and frequency. 

In Figure [7|-f8| it is clearly evident that the geometrical 
effects introduced by the Wilson depression are indeed sig¬ 
nificant for waves travelling towards the “umbra” {(j) < 50°), 
with the reduction in path length seemingly overriding the 
effects of the cooler plasma, resulting in similar (negative) 
travel-time shifts to those produced by the magnetic sunspot 
model. The combination of longer path length and cooler 
plasma results in positive St for waves travelling away from 
the spot centre, where in the magnetic sunspot model we 
observed negative Sr. A closer look at comparison plots in 
Figure|^also reveals that the travel time increase due to fast- 
to-Alfven conversion, typically seen around (j) ~ 80° — 100°, 
is absent in the thermal travel times. This indicates that the 
phase shifts produced by fast-to-Alfven mode conversion are 
indeed distinguishable from thermal/geometrical effects and 
have a distinct and significant effect on helioseismic travel 
time measurements in sunspots. 

The measured thermal St at 3 and 5 mHz appear to 
reach their peak a.t 9 ~ 32° in Figure However, we must 
remember that 9 in the thermal calculations is purely rep¬ 
resentative of the source position, not actual field inclina¬ 
tion from vertical, as magnetic effects are suppressed for 
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Figure 7. One-way phase travel time perturbations {St) derived 
from a model where only the thermal perturbations corresponding 
to the 1.5 kG sunspot with a Wilson depression of 400 km are 
present. Panels (a-b) represent a wave travel distance of A = 6.2 
Mm, (c-d) represent A = 8.7 Mm, and (e-f) represent A = 11.6 
Mm. Left column represents 3 mHz travel times and right column 
5 mHz. 


these calculations. Hence, this apparent dependance on 9 is 
a purely geometrical effect and is distinctly different from 
the “ramp effect” and fast-to-slow mode conversion-induced 
phase shifts we can observe in the magnetic St at 3 and 5 
mHz in Figure]^ 


6 DISCUSSION AND CONCLUSIONS 

As solar imaging hardware becomes increasingly sophisti¬ 
cated, the need for innovative diagnostic tools and precise 
modelling of wave propagation and transformation proper¬ 
ties in strong magnetic field regions of sunspots is made 
ever more apparent. Building on previous numerical stud¬ 
ies which employed simple plane-parallel atmospheres, we 
conducted a non-exhaustive parametric study of waves in 
model sunspot atmospheres in an attempt to further our 
understanding of the implications of MHD mode conversion 
on helioseismic measurements. 

By using time-distance heliosiesmology and a travel- 
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computations was also performed on the gSTAR national 
facility at Swinburne University of Technology. gSTAR is 
funded by Swinburne and the Australian Government’s Ed¬ 
ucation Investment Fund. Dr Shelyag is the recipient of an 
Australian Research Council’s Future Fellowship (project 
number FT120100057). 
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